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Abstract 

The paper (in French) exemplifies graphically a solution of the heat 
equation which is a 1-dimensional unfolding of an elliptic umbilic catas¬ 
trophe. The example is due to James Damon and adapts Thom-Mather’s 
singularity theory to multiscale models of scale-space analysis in image 
processing. 

Nous illustrons graphiquement un exemple, du a James Damon, de so¬ 
lution de I’equation de la chaleur qui est un deploiement a un parametre 
d’un ombilic elliptique. II fait partie d’une adaptation de la theorie des sin¬ 
gularites de Thom-Mather aux modeles multiechelle de I’analyse d’images 
(scale-space analysis). 


1 Introduction 

Pour evaluer I’heritage de I’oeuvre de Rene Thom, il faut je crois voir les choses a 
un double niveau. II y a d’abord la posterite directe de ses idees et puis il y a en- 
suite le developpement, dans d’autres contextes, d’idees de type thomien concer- 
nant I’application a des domaines originaux d’un certain nombre de formalismes 
de geometrie differentielle, de theorie des singularites, de theorie des bifurca¬ 
tions, de ruptures de symetrie, de phenomenes critiques, de phenomenes de tran¬ 
sitions de phases qui, jusque-la, ne concernaient que la physique. Rene Thom a 
ete I’un des principaux fondateurs d’un nouveau paradigme, le paradigme que 
j’appelle “morphodynamique” (theorie dynamique des formes), pour la biologie 
theorique et les sciences cognitives et ce dernier est devenu dominant dans de 
nombreux domaines. 

En ce qui concerne la posterite directe des idees de Rene Thom, il y a eu 
bien sur le travail de ses disciples. Meme si Rene Thom n’a pas constitue 
d’ecole a proprement parler, il a feconde la pensee d’un nombre considerable 
de personnes, aussi bien dans le domaine des mathematiques pures que dans 
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celui des applications et celui de la philosophie. Tons les chercheurs qu’il a 
inspires ont continue a developper certains aspects de son oeuvre meme si, sur 
le plan sociologique et mediatique, apres les grands debats des annees 70, la 
theorie des catastrophes congue an sens strict a connu un certain retrait. Mais 
ce retrait, d’ailleurs tout relatif, a ete, selon moi, tres largement compense par 
I’extraordinaire succes d’idees de type thomien. En particulier dans le domaine 
des sciences cognitives que je connais un pen, les idees introduites par Thom et 
Zeeman a la fin des annees 60 sont devenues des idees-force. 

L’un des centres d’interet majeurs de Rene Thom a ete d’elaborer une theorie 
des formes qui puisse etre compatible a la fois avec la genese objective de mor¬ 
phologies a I’interieur des substrats materiels physico-chimiques et biologiques 
(cf. ses travaux fondamentaux sur I’embryogenese) et aussi avec la perception 
subjective de morphologies a partir du traitement cognitif du flux optique man- 
ifestant ces substrats0 

L’un des premiers principes de la metaphysique thomienne de la nature 
est que les singularites constituent, en tant que discontinuites qualitatives, 
I’interface phenomenologique entre le monde objectif et la conscience percep¬ 
tive. Cela etait deja bien connu pour I’espace depuis I’esthetique transcendan- 
tale kantienne, I’espace operant a la fois comme cadre de I’objectivite physique 
et comme comme forme de la manifestation phenomenale, mais Rene Thom a 
ajoute a ce fondement de I’experience sensible I’intuition des discontinuites qual¬ 
itatives comme brisures de symetrie de I’homogeneite spatiale et cas generalises 
de transitions de phase (Thom [1972], [1980]) 11 

Dans mon hommage a Rene Thom de 1994 dans I’ouvrage Passions des 
formes edite par Michele Porte, j’insistais beaucoup sur I’actualite scientiflque 
des theses thomiennes dans les sciences cognitives, en particulier en ce qui con- 
cerne les structures de la vision. J’aimerais aujourd’hui dans cette courte note 
commenter un exemple allant non pas de la theorie des singularites vers la 
theorie de la vision mais de la theorie de la vision vers la theorie des singu¬ 
larites. 

En effet, comme pour toutes les idees seminales de Rene Thom, I’approfon- 
dissement du cote perceptif de la theorie morphologique a conduit a des problemes 
mathematiques originaux. Pour comprendre le role des singularites dans la per¬ 
ception des formes, il a fallu faire la synthese entre la classification de Thom- 
Mather et le point de vue naturel de I’analyse geometrique des images en vision 
computationnelle, qui est celui des modeles dits multi-echelle. La synthese a 
ete effectuee par James Damon dans une serie d’articles et en particulier dans 
le papier “Local Morse theory for solutions to the heat equation and Gaussian 
blurring” (Damon, 1995) sur lequel nous nous appuyerons ici. 

La classification des modeles multiechelle de singularites fait apparaitre cer¬ 
tains phenomenes nouveaux dont le plus surprenant est celui de la possibilite 

^Le lecteur interesse par les differents aspects de la theorie de la forme chez Rene Thom 
pourra consulter notre article “Forme” dans VEncycloprEdia Universalis (Petitot [1989]). 

^Pour une introduction aux multiples facettes de la philosophie de Rene Thom, on pourra 
se referer aux actes du Colloque d’hommage tenu a Cerisy en 1982 ainsi qu’a leur presentation 
synthetique dans Physique du Sens (Petitot [1992]). 
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de creation de points critiques lorsque I’echelle augmente, evenement qui sem- 
ble violer le principe du maximum auquel satisfont ces modeles. D’oii I’interet 
de preciser ce phenomene curieux. Nous le ferons ici au moyen de graphiques 
calcules avec Mathematica™. 

Ces petits calculs ont ete efFectues il y a quelques annees lors d’un colloque 
a, la Fondation des Trebles organise avec Bernard Teissier, Jean-Michel Morel 
et David Mumford dans le cadre du Trimestre special du Centre Emile Borel 
“Questions mathematiques en traitement du signal et de I’image” (septembre- 
decembre 1998). Ils devaient servir de base a des discussions avec Rene Thom 
mais la vie en a decide autrement. 


2 Le concept de “Scale space analysis” 

L’idee d’analyse multiechelle, de “scale space analysis” ou d’algorithmes “pyra- 
midaux”, domine les theories de I’analyse geometrique des images. Elle remonte 
a Witkin (1983) et a Koenderink (1984, 1986) et part de la remarque suivante. 
Pour etre morphologiquement correcte, une analyse des images doit s’effectuer 
en termes de geometrie differentielle. Le probleme est que les outils de cette 
derniere ne sont pas directement applicables au signal en tant que tel, qui est 
trop bruite pour etre differentiable. Pour que les images puissent acquerir le 
statut d’observables geometriquement analysables par detection d’invariants, il 
faut par consequent definir au prealable une echelle, c.a.d. fixer un niveau de 
regularisation du signal bruite. Jan Koenderink I’a souvent souligne : le scaling 
est essentiel. C’est ce que fait d’ailleurs la vision naturelle puisque les neurones 
visuels primaires possedent un champ recepteur qui est un petit domaine de 
photorecepteurs retiniens et operent sur le signal optique comme des filtres par 
convolution avec leur profil recepteur (leur fonction de transfert) qui est une 
fonction definie sur leur champ recepteur. 

Mais comment une analyse multiechelle peut-elle deboucher sur une veritable 
analyse morphogenetique d’une image 2D definie comme une fonction (ou, mieux, 
une distribution) I{x,y) sur la fenetre retinienne R de coordonnees x et y? 
L’idee directrice est de plonger I’image dans une famille Is{x, y) parametree par 
une echelle s de fagon a ce que : 

(i) lo = I, 

(ii) Il soit une image indifferenciee, et 

(iii) Lorsque I’echelle s croit, revolution de Iq a Ji “simplifie” strictement 
I’image. Cette contrainte dite de “causalite” interdit I’apparition ex ni- 
hilo de nouveaux details lorsque I’echelle croit. 

L’evolution avec I’echelle s des lignes de niveau de Is, c’est-a-dire la suite 
d’evenements de bifurcation qu’elles subissent en se simplifiant progressivement, 
fournit une methode puissante pour analyser la structure morphologique de 
I’image et sa decomposition en elements constituants. On montre que, sous 
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des contraintes generales de linearite, d’invariance par translation, d’isotropie 
et d’invariance d’echelle, la fagon la plus simple d’obtenir un tel resultat est de 
prendre pour Ig une solution de Vequation de dijfusion typique qu’est I’equation 
de la chaleur dlg/ds = Alg. Dans la mesure on le noyau de la chaleur est 
gaussien, on est ainsi conduit a I’idee d’un lissage gaussien multiectielle de 
I’image (Gaussian blurring)!^ 

Mais ce qui domine phenomenologiquement une image sont les singularites. 
Si Ton veut adapter a I’analyse d’images la theorie de Thom-Mather des sin¬ 
gularites generiques et des deploiements universels, il faut par consequent la 
rendre multiechelle. 


3 Theorie multiechelle des singularites 

James Damon a montre comment on pouvait transformer la theorie de Morse- 
Whitney-Thom-Mather-Arnold dans le cas de I’equation de la chaleur. La 
difhculte principale est que les formes normales de Morse ne satisfont pas a 
cette EDP. Les methodes doivent done etre transposees des espaces de germes 
d’applications C°° aux espaces de germes de solutions de I’EDP, or ces espaces 
n’ont pas les “bonnes” proprietes algebriques qui font marcher la theorie et per- 
mettent d’appliquer des theoremes de transversalite de Thom pour obtenir des 
resultats de genericite0 

James Damon a done du d’abord redefinir le concept d’equivalence pour des 
germes (7°° d’applications f{x, s) 

f : (R"+i,0) —(R,0) 

(x,s) I—;> z = f{x,s) 

oil X = {Xi, . . . ,Xn) ■ 

• f et g sont iL-equivalentes s’il existe un germe de diffeomorphisme ip : , O) —>■ 

(K"+^,0)de la forme ip{x,s) = s), (p 2 (s)) avec > 0 et un germe 

c : (R, 0) —>■ (R, 0) tels que 


9 {x,s) = fo(p{x,s) -bc(s). 

Autrement dit, la difference entre g et f peut etre resorbee 

(i) par des translations au but c(s) dependant de s, 

(ii) par un changement de coordonnee (preservant I’orientation) ip 2 sur I’axe s 

de la source, 

®I1 faut insister sur le fait que s est ici un parametre d’echelle et non pas un parametre 
temporel comme e’est le cas d’habitude. L’evolution s’effectue dans un espace-echelle et non 
pas dans un espace-temps. 

“^Pour une premiere introduction a la theorie des singularites de Thom-Mather on pourra 
consulter en particulier le seminaire Bourbaki Chenciner [1973], Particle de VEncyclopcEdia 
Universalis Chenciner [1980] et le Graduate Text Golubitsky&Guillemin [1973]. Les six arti¬ 
cles de base de John Mather sont inclus dans la bibliographie. 
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(iii) par des changements de coordonnees dependant de s, du sous- 

espace R" de la source. 

• De meme f&tg sont /S'-equivalentes {IS signifie “Intensity Sensitive”) si c 
est constantly et s’il existe de plus un germe de diffeomorphisme tp : (M^, O) —>■ 
(R^,0) de la forme ip{z,s) = {6(z,s),s) avec f|(0, 0) > 0 et 0(0, s) = 0 pour 
tout s tels que 


g{x,s) = 9 o f o (p{x,s) +c. 

Cela signifie que Ton peut remplacer c(s) par des changements de coordonnees 
appropries au but (preservant I’orientation) 9{», s) dependant de s. 

A ces notions d’equivalence sont naturellement associees des notions de sta- 
bilite par deformation. On considere des families fw{x,s) de fonctions f{x,s) 
parametrees par des parametres li; G 

U: (R"+i+9,0) ^ (R,0) 

{x,s,w) I—^ z = f{x,s,w) = fwix,s) 

et deformant f{x, s, 0) = fo{x, s). 

Definition. fo{x,s) est iJ-stable (resp. IS'-stable) si fQ{x,s) est son pro- 
pre deploiement universel, autrement dit si toute deformation assez petite est 
triviale pour la 7J-equivalence (resp. la /S'-equivalence). 

En adaptant les demonstrations de Rene Thom et de John Mather, James 
Damon a trouve les formes normales des singularites stables et des deploiements 
universels de bas degre. Par exemple les equivalents iJ-stables des points 
critiques quadratiques (points critiques non degeneres stables ) pour la H- 
equivalence ont pour forme normale : 

i—n 

±^xl± (2n)s 
et 


i—n i—n 

aiX^ avec ^ 0 pour tout i et at = 0 

i=l i=l 

dont on verifie trivialement qu’ils satisfont I’equation de la chaleur. Mais on 
obtient aussi des singularites iJ-stables cubiques de forme 

xl + 6sxi + Q{X 2 , ..., Xn, s) 
et 


xl — 6sa:i H— 6xixl + Q{x 2 ,..., Xn, s) 

= 0 si ^( 0 , 0 ) = 0 , mais on peut avoir a considerer des cas c 7 ^ 0 . 
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oil Q{x2, ... ,x„,s) est une singularite quadratique en {x2, ...,x„). Cela est du 
au fait que la variable d’echelle s est de poids 2 relativement aux variables Xi 
qui sont de poids 1 et qu’un terme apparemment quadratique comme sxi est 
done en fait cubique. 

Pour la /S'-equivalence, on obtient la meme premiere forme de singularites 
quadratiques stables 


i=n 

i=l 

Mais la seconde forme n’est plus la meme et devient 

i—n / i=n \ i—n 

aiX^ H- 2 j ai ] s avec ^ 0 pour tout i et ^ 0 

\i^l ) 

i—n 

Pour ^ = 0 on obtient une nouvelle forme plus compliquee qui, dans le cas 

i—1 

n = 2, s’ecrit : 



0. On obtient egalement deux formes cu- 


s avec Oi ^ 0 pour tout i et 


avec Ui ^ 0 pour tout i et ai 

i—1 

biques /S'-stables : 


cf + 6sxi + E + 2 E 


et 


i—n / i—n \ i—n 

Xj—6sxiH—6X1X2+E aiXi +2 ( E 1 ® ai 0 pour tout i et E«*^o 

i—2 \i—2 / i—2 

C’est cette derniere forme qui est surprenante car elle conduit a des creations 
de points critiques, ce qui pourrait sembler interdit par le principe de causalite 
expose plus haut. 


4 Analyse graphique de la creation de points cri¬ 
tiques par diffusion 

Nous allons analyser graphiquement la creation de points critiques dans le cas 
n = 2. Nous etudions done la solution de I’equation de la chaleur : 
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Figure 1: Les points critiques de f{x,y,s) se regroupent en 2 branches. 


f{x, y, s) = x^ — 6xy'^ + y'^ — 6sa; + 2s 

consideree comme une famille a un parametre (I’echelle s) de fonctions de {x, y). 


4.1 Forme des points critiques 

Les points critiques sont donnes par 


— 6y^ — 6s = 0 
II = -l2xy + 2y = Q 


dont les 4 solutions sont 


(l-72s)^/^ 


xi± - \ , yi± - ^ 

X2± = ±\/^, J/2± = 0 


Ils constituent deux branches paraboliques parametrisees par s, pci(s) et pc 2 (s). 
Pour s < 0 (ce qui est interdit pour une echelle), pc 2 {s) est imaginaire. Pour 
s > 1/72 c’est la branche pci qui devient imaginaire, pc 2 donnant deux points 
critiques opposes. Pour 0 < s < 1/72, il y a 2 paires de points critiques (indexees 
par 1± et 2±) : la branche pc 2 apparait (creation de deux points critiques 
opposes) et reste en competition avec pci tant que s < 1/72. Cf. Figure [T] 

Le Hessien H de f etant donne par 


H = 


dxidxi J 


sa valeur aux points critiques est 


6x —12y \ 

-12y 2- 12x )' 




1 _ Tv'2(1-72s) 

Tv' 2(1 - 72s) 0 


et H 2 ± 


±dVTs 0 

0 2 (1 T 6^) 


7 








et ses valeurs propres (dans I’ordre de i^l+, Hi-, H 2 +, H 2 -) sont 


f 1 - 3V1 - 64s l + SVl-eds] f 1 + 3V1 - 64s 1 - 3V1 - 64s 
I 2 ’ 2 J ’I 2 ’ 2 

|2 - 12v^, 6v^|, |2 + 12v^, -6v^|. 

Quant aux vecteurs propres, on trouve, toujours pour le meme ordre. 



j-1 + 3V1-64S 
\ 2i/2(l - 72s) ’ 

[ 1 - 3vT^64i / 

\ 2^2(1-725)’ j’\ 


-1 - 3V1 - 64s 
2i/2(l - 72s) ’ 

1 + 3V1 - 64s 



2a/ 2(1 - 72s) 

{{ 0 , 1 },{ 1 , 0 }},{{ 0 , 1 },{ 1 , 0 }}. 

Comme les deux points critiques de la branche pci deviennent imaginaires 
conjugues pour s > 1/72, ce sera a fortiori le cas pour s > 1/64. Les change- 
ments de signe des valeurs propres des points critiques reels s’effectuent pour 
1 — iy/l — 64s = 0 et 1 — = 0, ce qui dans les deux cas donne s = 

D’ou les tableaux suivants ou A et /x sont les deux valeurs propres. 

(i) Pour 0 < s < 1/72 : 


(1 


pci+ = 

(l-72s)i/= \ 
6V2 ) 

A < 0 

PCl- = 
(1 (l-72s) 

V6’ 6\/2 

A > 0 

1 / 2 ^ 

PC2+ = 
(v^,0) 

A > 0 

PC2- = 
-(v^,0) 
A > 0 

/X > 0 

col 

^ < 0 
col 


^ > 0 
sommet 

^ < 0 
col 

1/72 : 






PC2+ = 

(v^, 0) 

A < 0 

PC2- = 

(v^,0) 

A > 0 



/X > 0 
col 

/X < 0 
col 




4.2 Representations des surfaces 

Aux points critiques les valeurs de f{x,y,s) sont zi± = s + 1/216, Z 2 ± = 
2s (1 =F 2-\/2s). Pour bien visualiser la famille de fonctions f{x, y, s) nous aliens 
representer de plusieurs manieres les graphes de leurs sections a s = cste. 



La planche 1 represente 6 sections pour s variant de —(1/3)1/72 a (4/3)1/72 
par pas de (1/3)1/72. Le point de vue par le dessus permet de bien visualiser 
la position des points critiques. 




Planche 1 

La planche 2 represente les memes sections mais dans un point de vue plus 
frontal permettant de mieux voir le relief de /. 
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Planche 2 

Quant a la planche 3 elle montre les lignes de niveaux de / ainsi que son 
champ de gradient V/ par rapport aux variables spatiales x et y. 
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La figure [Hzoome sur la situation s = (2/3)1/72 proche de la valeur critique 
s = 1/72. 

On voit bien la suite des bifurcations. Au depart (s < 0) il n’y a que les 
deux cols pci± symetriques par rapport a I’axe des x et d’abscisse fixe 1/6. Ils 
se rapprochent de I’axe des x comme i/l — 72s. Pour s = 0 apparaissent sur 
I’axe des x un sommet pc 2 + d’abcisse et un colpc 2 - d’abscisse —Pour 
s = 1/72, pci+ et pci- fusionnent sur I’axe des a; en a: = 1/6 au moment oil 
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Figure 2: Zoom sur la valeur de I’echelle s = (2/3)1/72 proche de la valeur 
critique s = 

pc 2 + est precisement en x = \/^ = 1/2/72 = 1/6. La fusion de ces trois points 
critiques transforme pc 2 + en un col qui s’eloigne sur I’axe des x en -/Ss- 

On voit clairement sur ces images I’apparition temporaire du maximum. Si 
le principe de causalite n’est pas viole c’est a cause de la fagon dont les surfaces 
se positionnent les unes par rapport aux autres. La nappe montante reste a peu 
pres fixe alors que la nappe descendante s’abaisse comme le montre la figure [3] 
qui represente la section le long de I’axe de symetrie y = 0. 



Figure 3: Section des surfaces f{x, y, s) le long de I’axe de symetrie y = 0. 
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Figure 4: La section du deploiement universel de I’ombilic elliptique gw,u,v = 
— 6xy'^ + wx'^ + ux + vy pour w = 1/2. 

4.3 Plongement dans le deploiement de I’ombilic ellip¬ 
tique 

La famille f{x,y,s) = x^ — Gxy"^ + y^ — dsx + 2s est une deformation de la 
singularite ombilic elliptique go{x,y) = x^ — 6xy^ et il est done pertinent de 
la positionner par rapport an deploiement nniversel standard de ce dernier, 
gw,u,v = x^ — 6xy'^ + wx^ +ux + vy. En fait e’est ce qui se passe an voisinage dn 
point critique degenere triple (1/6,0) pour s = 1/72 qui est le plus interessant. 
On a 


/(l + i,j,,.)=i=>-fiV + y+(^-6s)x + »+2L 

ce qui, pour s = 1/72, donne la singularite (centre organisateur) 

./I 1\ ^ 2 1 

Le deploiement universel gw,u,v = x^ — 6xy^ + wx^ +ux + vy a pour discrimi¬ 
nant dans le {u, v, w) une surface conique sur une hypocycloi'de a 3 rebrousse- 
ments. Sa section plane w = 1I2 est une telle hypocycloi’de H d’equation 

(-1 -b u)u^ + (486 - 648m -b lUu^)v‘^ + 5184n^ = 0 
Elle est representee a la figure IH 

Pour se faire une idee de la structure de I’ombilic elliptique on pent considerer 
dans le {u,v,z) le graphe des valeurs critiques de gi/ 2 ,u,vix,y). Les figures 
[S]et[n] en donnent deux points de vue. 

Pour visualiser le detail de la structure pres de I’origine, on a represente a 
la planche 4 trois sections au-dessus de trois fenetres du plan de controle (it, v). 
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Figure 5: Un point de vue sur le graphe des valeurs critiques de Tombilic ellip- 
tique gi/ 2 ,u,v{x,y). 



Figure 6: Un autre point de vue sur le graphe des valeurs critiques de rombilic 
elliptique gi/ 2 ,u,v{x,y). 


14 





Figure 7: Le positionnement de la solution f{x, y, s) de I’equation de la chaleur 
dans le deploiement de I’onibilic elliptique. 






On voit que la solution stable de I’equation de la chaleur f{x, y, s) corre¬ 
spond a un deploiement a un parametre hautement non generique de I’ombilic 
elliptique. Dans le (m, v, c) le discriminant du deploiement gw,u,v = x^ — 
Gxy"^ + + ux + vy + c est le cylindre de base H. La solution stable f{x, y, s) 

correspond a la droite u = (^^ — 6s), c = s -I- ^ au-dessus de I’axe u qui rentre 
dans H pour s = 0 en u = 12 a la hauteur c = 1/216 et en ressort par le cusp 
u = 0, c = 1/54 pour s = 1/72 (cf. Figure[7]). 
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5 Equations de diffusion non lineaires et traite- 
ment d’images 


Nous avons vu comment la theorie des singularites de Whitney-Thom-Mather- 
Arnold pouvait etre adaptee aux solutions d’une EDP parabolique de diffusion 
comme I’equation de la chaleur. Mais il ne s’agit la que du debut de Thistoire. 
En effet les methodes de diffusion simple (homogene et isotrope) font probleme 
car le lissage (Gaussian blurring) qu’elles effectuent est indifferent a la geometrie 
de I’image. D’oii I’idee qui s’est progressivement imposee au cours des annees 
80 que, pour pouvoir effectuer une bonne analyse morphologique des images, il 
faut concilier deux exigences apparemment contradictoires : 

(i) regulariser le signal de fagon multiechelle par diffusion; 

(ii) preserver la geometrie de I’image, c’est-a-dire les discontinuites d’elements 

differentiels possedant une signification geometrique intrinseque. 


Pour ce faire, il faut adapter I’equation de diffusion a la preservation de ces 
elements differentiels construits a partir des jets successifs de la fonction /^(x, y). 

Le cas le plus simple est celui des bords qui delimitent les domaines ho¬ 
mogenes d’une image et sont essentiels a la definition de ses constituants. Un 
bord est idealement une discontinuite du gradient V/s(x,y) de Ig. Pour qu’une 
diffusion preserve le caractere discontinu des bords tout en les simplifiant pro¬ 
gressivement, il faut qu’elle soit anisotrope et, plus precisement, inhibee dans 
la direction du gradient. L’EDP la plus simple possedant cette propriete est 
I’equation de diffusion non lineaire : 


dis 

ds 


|V/s| div 


\VIs 


= AA - 


g(V/„VA) 

|VA|" 


oil H est le Hessien de Ig. Elle est uniformement parabolique le long des combes 
de niveau de Ig mais totalement degeneree dans la direction du gradient. Elle 
fait evoluer les lignes de niveau — et done en particulier les bords — comme 
des fronts avec une vitesse normale egale a leur courbure. 

Jean-Michel Morel et ses collegues (Alvarez et al. 1992) ont introduit un 
lissage gaussien supplementaire dans ces EDP de fagon a controler la vitesse de 
diffusion en la couplant a la geometrie de I’image. L’equation de base devient 
alors : 


A = 9(|G.V/.|)|V/.|div(^) 

oil G est un noyau gaussien et g{x) une fonction decroissante telle que g{x) —> 0 
quand x —>■ oo. Ils ont egalement montre que si I’on impose une contrainte 
d’invariance affine (au lieu de I’invariance euclidienne), I’EDP type devient : 


dig 

ds 


|VA| div 


f ^Is \ 
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On pent generalise! ces equations de fagon a ameliorer I’analyse morphologi- 
que de I’image. Par example les “cretes” d’une forme — qui fournissent une 
“squeletisation” essentielle a son analyse morphologique, ce que Rene Thom 
avait souvent souligne a la suite de Harry Blum en utilisant le concept de “cut- 
locus” ) — peuvent etre extraites du signal en tant que des discontinuites de la 
direction du gradient qui sont preservees de la diffusion. 

En considerant les lignes de niveau des fonctions Is{x,y) on obtient des 
evolutions de courbes planes Cg fermees qui se propagent comma des fronts 
conformement a une loi du type ^ = F{K)n ou p est un point de Cs, n la 
normale (externe) en p a Cg et k la courbure de Cg en p . Les cas les plus 
etudies sont : 

(i) celui de la propagation a vitesse constante : || = ±n, ^ (modeles 

de propagation d’ondes de type “grassfire” engendrant le cut-locus de la 
courbe), et 

(ii) celui de la propagation a une vitesse proportionnelle a la courbure (r est 

I’abscisse curviligne de la courbe) : || = —Atn, ^ 

(iii) Stanley Osher et James Sethian (1988) ont egalement etudie les cas in- 

termediaires ou F{k) = 1 — ek. La courbure k satisfait alors une equation 
de reaction-diffusion du type : ^ -I- ek^ — k^. 

Sous les titres de “curve shortening”, “flow by curvature” ou “heat flow 
on isometric immersions”, le cas (ii) a ete particulierement investigue par des 
geometres comme M. Gage, R. Hamilton (1986), M. Grayson (1987), S. Osher, 
J. Sethian (1988), J. Sethian (1990), L.G. Evans et J. Spruck (1991)@ Si jg : 

^ est I’immersion isometrique definissant Cg, comme on a Ajg = —ku, 
I’equation de diffusion ^ = —nn est en fait I’equation de la chaleur (pour les 
immersions) ^ = Ajg. Dans I’espace fonctionnel J des immersions j : ^ 

cette equation definit le champ de gradient de la fonction de Morse donnant 
la longueur de la courbe image C = j{S^). 

Un theoreme fundamental de “curve shortening” du a M. Gage et R. Hamil¬ 
ton (1985) et generalise par Matthew Grayson (1987) dit que si Co est une 
courbe plongee dans le plan (meme tres sinueuse), alors I’equation de la chaleur 
la contracte sur un point, Cg devenant asymptotiquement un cercle (conver¬ 
gence pour la norme C°°). En particulier Cg devient convexe avant de pouvoir 
developper des singularites. 

®En fait, la theorie vient de Richard Hamilton qui cherchait a r&oudre des problemes de 
Relativite generate. En utilisant I’equation de la chaleur, il a montre que si X est une variete 
riemannienne compacte de dimension 3 avec courbure de Ricci Rij > 0, alors X admet une 
metrique riemannienne a courbure > 0 constante. Or, ces dernieres sont classifiees. II cherchait 
egalement a engendrer des geodesiques fermees a partir de courbes fermees quelconques. (Cf. 
le seminaire Bourbaki Bourguignon [1985]). 
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6 Conclusion 


C’est en fait toute la conception morphodynamique thomienne sur le role des 
singularites, des bifurcations et des ruptures de symetries dans la perception 
qui se trouve confirmee par un nombre croissant de travaux contemporains en 
vision computationnelle et en neurosciences de la vision. Citons pour conclure 
trois examples : 

(i) Celui des processus de segmentation des images en region liomogenes delimitees 

par des discontinuites qualitatives. Ce probleme central etait a la base des 
modeles morphologiques de Thom des le debut avec la difference entre les 
points reguliers d’un substrat ou toutes les qualites phenomenologiques 
(couleur, texture, etc.) sont localement homogenes et les points singuliers 
ou au moins une qualite phenomenologique subit une discontinuite qual¬ 
itative. II a ete formule comma un (difficile) probleme variationnel de 
type “free boundary problem” par David Mumford dans les annees 80 et 
a donne lieu a tout un ensemble de travaux, en particulier autour d’Ennio 
De Giorgi a Pise (Luigi Ambrosio, Gianni Dal Maso). La conjecture de 
Mumford sur la nature des points singuliers des solutions optimales a 
presque ete demontree par Alexis Bonnet et Guy David. 

(ii) La squeletisation des formes planes au moyen de leur cut locus, c’est- 

a-dire des singularites d’un processus de propagation a la Huyghens a 
partir du bord. Cette idee profondement morphologique a ete introduite 
dans les annees 60 par le specialiste de la vision Harry Blum [1973] et 
a toujours ete tres fortement defendue par Rene Thom. L’avantage du 
squelette d’une forme est d’etre un graphe de dimension 1 constitue d’arcs 
se rejoignant generiquement en des points triples ou s’arretant en des 
points singuliers. II fournit done automatiquement une decomposition 
qualitative de la forme en constituants qui sont des cylindres generalises 
(cf. Marr [1982]). L’homeomorphismes des squelettes definit une relation 
d’equivalence qui est plus forte que les diffeomorphismes et plus faible que 
les isometries tout en etant tres differente de I’equivalence conforme. Pen¬ 
dant longtemps cette methode etait peu prisee car instable et trop sensible 
aux deformations du bord, toute petite asperite du bord engendrant une 
branche supplementaire du squelette. Mais le developpement de methodes 
de pruning multiechelle a rendu la methode stable et elle est aujourd’hui 
tres utilisee. Un des meilleurs specialistes en est Benjamin Kimia. 

(iii) La reconstruction des surfaces plongees dans a partir de la famille de 
leurs contours apparents. La encore Rene Thom a toujours insiste sur 
I’importance du probleme. Si S est une surface dans a chaque point 
de vue p (plan de projection et direction de projection, les p torment une 
grassmannienne V) est associe un contour apparent Cp. Le type qualitatif 
de Cp stratifie V a travers une hypersurface H et le probleme est de 
comprendre I’equivalence entre la geometrie de S et la geometrie de la 
stratification definie par H. Par exemple si Ton connait des equations 
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algebriques pour S peut-on en deduire des equations pour H? Le probleme 
est difficile. Or le systeme visuel est capable d’anticiper avec une precision 
remarquable les variations de Cp lorsque Ton tourne autour de S. On ne 
connait pas les algorithmes neuralement implementes sous-jacents a cette 
performance perceptive 0 

Tons ces travaux convergents montrent que le reve de Rene Thom de developper 
une “physique phenomenologique” et une “semiophysique” de la perception est 
en train de se realiser. II en va de meme dans le domaine de la morphogenese 
biologique ou Rene Thom estimait avoir introduit ses idees les plus fondamen- 
tales. La aussi, comme pour les sciences cognitives, on pent dire que I’utilisation 
de modeles dynamiques et de la theorie des bifurcations engendrant des discon- 
tinuites qualitatives dans des substrats physiques, s’est beaucoup developpee 
parallelement a d’autres modeles de theorie des patterns comme les modeles de 
reaction-diffusion de type Turing. Ce qui a peut-etre le plus change d’avec la 
fagon dont Rene Thom lui-meme voyait les choses, est qu’il semble desormais 
possible d’etablir une convergence avec les resultats de la genomique. L’on com¬ 
mence a comprendre le controle genetique du developpement et des reactions 
morphogenetiques qui interessaient Rene Thom. 

Rene Thom etait non seulement un puissant genie mathematique et un 
maitre humainement extraordinaire, mais aussi Pun des esprits philosophiques 
les plus inspires du XXe siecle. Ce fut un privilege inoubliable que d’avoir pu 
collaborer etroitement et continument avec lui. Au-dela des aleas humains trop 
humains de la sociologie des milieux academiques, ses idees ont d’ores et deja 
triomphe dans de nombreux domaines et, j’en suis convaincu, continueront a le 
faire, ainsi que toutes les grandes idees appartenant au patrimoine de la pensee. 
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